home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48hor1 / planets.doc < prev    next >
Text File  |  1995-03-31  |  5KB  |  114 lines
  1. (Comp.sys.handhelds) 
  2. Item: 3528 by TDSTRONG%MTUS5.cts.mtu.edu at VM1.NoDak.EDU 
  3. Author: [Tim Strong] 
  4.   Subj: Planetary postions HP48SX 
  5.   Date: Thu Jun 20 1991 
  6.  
  7. Hello, 
  8.  
  9. Here is the program I promised about a week ago.  It calculates the postions 
  10. of the planets in the sky. 
  11.  
  12. You don't have to know anything about julian days or anything like that; 
  13. the planet finder invokes them internally.  However the functions are 
  14. available. 
  15.  
  16. Simply load the directory into your calculator and press CST.  (Incidentally 
  17. make sure your machine is in DEG mode or your answers will come out in 
  18. radians). After pressing CST you will get the following options: 
  19.  
  20. PLAN            This key runs the planet finder.  I have tried to make it as 
  21.                 user friendly as possible since after a few weeks I'd forget 
  22.                 how to use it if it weren't. 
  23.  
  24. CONC            This key invokes a menu driven astronomical coordinate 
  25.                 conversion program.  I will describe its use in a few minutes. 
  26.  
  27. DMS->   This key takes two coordinates that are on the stack in 
  28.                 HH.MMSS or DDD.MMSS format and converts them to decimal. 
  29.                 Basically it invokes the built-in HMS-> twice. 
  30.  
  31. ->DMS   This key does the reverse of the above but also tags the output 
  32.                 so you know it's in HH.MMSS form. 
  33.  
  34. JULDAY  This key finds the julian day for any date.  Its input is month 
  35.                 on stack level three, day on level two, and year on level 1. 
  36.  
  37. deltaD90        This key finds the number of days since the epoch of 1990, 
  38.                 which is the epoch the internal coordinates or values for the 
  39.                 planets orbit elements are set for. 
  40.  
  41. epoch   This variable contains the current epoch in which you want the 
  42.                 coordinate to be figured.  If your star maps are labled epoch 
  43.                 2000 you put 2000 in here. 
  44.  
  45. phi             This variable contains your local latitude in decimal form. 
  46.                 Set it for your viewing location. 
  47.  
  48. epsilon   This variable will find the obliquity of the eliptic for 
  49.                 for the current epoch.  [Say THAT five times fast! -jkh] 
  50.  
  51. First try PLANETS.  A menu across the bottom comes up showing all the planets 
  52. and you are instructed to "ENTER DATE M DD YYY and select planet".  Simply do 
  53. that (for example if I want to know all about Jupiter tonight I enter 6 20 1991 
  54. and then press JUPIT).  After ten seconds I get the following display (the 
  55. Greek letters are spelled out): 
  56.  
  57. alpha= 9.0000 
  58. delta= 17.4839 
  59. rho= 6.0239 
  60. tau= .5006         Jupiter 
  61. theta= 32.6598 
  62. F= .9956 
  63. m= -1.8574 
  64.  
  65. Alpha is the right ascension; it reads as 9hr 00min 00sec.  Delta is the 
  66. declination which reads as 17deg 48min 39sec.  Rho is the number of 
  67. astronomical seconds.  Theta is the angular diameter of the planet (how large 
  68. it will appear in a telescope) in arcminutes.  F is the current phase of the 
  69. planet, as percent lighted.  M is the visible magnitude of the planet. 
  70.  
  71. Now the numbers for right ascension and declination should be accurate to 
  72. within a minute or two (Mercury may be worse).  Also the magnitude is 
  73. probably a bit optimistic. 
  74.  
  75. Next press CONCO.  If there are less than three items on the stack you will see 
  76. a menu of conversions and the words "select conversion:".  Put the coordinates 
  77. you want to convert on the the stack in decimal notation (use ->DMS if you want 
  78. to) and then press the appropriate conversion.  If there are three items on 
  79. the stack the program assumes it has been called by another program.  In this 
  80. case the following arguments should be on the stack: 
  81.  
  82. 3: first coordinate 
  83. 2: second coordinate 
  84. 1: number of conversion 
  85.  
  86. The number of the conversion is simply the number of the conversion you want 
  87. from the menu numbered left to right. 
  88.  
  89. Some facts about the program. 
  90.  
  91.    1)  Please feel free to modify it in any way.  If you find a method that 
  92.        makes more sense, is shorter, faster, or easier to understand, or 
  93.           adds features please post it. 
  94.  
  95.    2)  The program is quite long 4679.5 bytes.  I did that to make it 
  96.        user friendly and because it's my first large program (it's not 
  97.        particularly efficient). 
  98.  
  99. The following reference should be cited as I used it heavily in order to 
  100. find an algorithm to find the values: 
  101.  
  102. Practical Astronomy with your Calculator--Peter Duffet Smith  3rd. edition 
  103. Published by the Cambridge University Press 1988. 
  104.  
  105. This book is excellent.  It even shows how you can use you calculator to 
  106. predict and picture a solar eclispe.  All the routines shown are algorithms 
  107. that you work out by hand, writing down intermediate results, but they can be 
  108. programmed with some work (kinda like I did 8-)) 
  109.  
  110. ------------------------------------------------------------------------- 
  111.    __ 
  112.   I__)  _   _I  _   _   Tim Strong <TDSTRONG@MTUS5.cts.mtu.edu> 
  113.   I  \ (_I (_I (_I I    Michigan Tech.    Houghton, Michigan, U.S.A. 
  114.